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Abstract. - Stunningly large exciton transfer rates in the light harvesting complex of photo- 
synthesis, together with recent experimental 2D spectroscopic data, have spurred a vivid debate 
on the possible quantum origin of such efficiency. Here we show that configurations of a random 
molecular network that optimize constructive quantum interference from input to output site yield 
systematically shorter transfer times than classical transport induced by ambient dephasing noise. 



Introduction. — A growing amount of experimental 
data [lj|3j suggest that quantum coherence may be at the 
origin of the stunning efficiency of exciton transport in 
photosynthetic light harvesting, even at ambient temper- 
atures and in a doubtlessly very noisy environment. The 
simplest natural structure that exhibits these surprising 
properties is the Fenna-Matthews-Olson (FMO) light har- 
vesting complex of green sulfur bacteria, which consists of 
seven (or eight [i] ) bacteriochlorophyll molecules arranged 
in a disordered network [5j[6j . Exciton transfer is here me- 
diated by dipole-dipole coupling between these different 
molecular sites, and is associated with a de-excitation of 
the donor site from some excited to the ground state, and 
an excitation of the acceptor site from the ground to some 
excited state. The molecular network that mediates the 
exciton transfer, from the antenna complex to the reaction 
center, where the excitation fuels the organism's chem- 
istry, is embedded into a complicated protein structure, 
which seems to provide some structural stiffness, and also 
defines a nontrivial spectral structure of the environment, 
which preserves the coherence on the FMO network itself, 
on considerably longer time scales than to be expected 
for a white noise environment [?]. Since the prevalence 
of the coherence effects in widely open systems in contact 
with high temperature environments challenges our tradi- 
tional understanding of what seemed the restricted realm 
of quantum mechanics, these highly specialized biological 
functional units let us think anew how to control efficient 
transport in disordered and noisy systems, possibly by ex- 



ploiting fundamental quantum features. 

The available experimental evidence, still not always 
fully consistent [l][2] and vividly debated, is mirrored by a 
large variety of theoretical approaches, which distinguish 
themselves in terms of the applied methodology as well 
as of the level of faithfulness of the modeling of the actual 
biomolecular object under scrutiny — from advanced quan- 
tum molecular dynamics [4|, over quantum simulations 



8|-|10 based on some effective, seven site Hamiltonian 
Lindblad equations [12| or various non-Markovian ap- 
proaches derived from open system theory 13 -17 , to di- 



verse quantum optical models 18 ■ 20 and abstract statis- 



tical treatments 21 22 . In essence, there are four princi- 
pal lines of argument to explain efficient exciton transport, 
by (i) noise-assisted, (ii) non-Markovian or (iii) driving, 
and (iv) multi-path interference induced exciton transfer. 
Given the available experimental data and the hitherto 
limited characterization and control of the precise micro- 
scopic Hamiltonian and environment coupling agents that 
generate the experimentally observed phenomena, it re- 
mains an open question which of these suggested mech- 
anisms were actually used by evolution to optimize the 
FMO's functionality. Since any improvement that pro- 
vides an evolutionary advantage will be implemented, it 
is even not unlikely that all of them are used at some 
level. However, in the light of the debate about a pos- 
sible quantum enhancement of transport in the photosyn- 
thetic complex, it is highly relevant to understand the spe- 
cific role of quantum coherence for these different mecha- 
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nisms, and to compare the achievable transfer efficiencies 
they aUow for. It is the present contribution's purpose 
to provide such comparison for noise-enhanced and multi- 
path-interference transport scenarios. Large scale statisti- 
cal sampling will allow us to show that multi-path quan- 
tum interference always leads to better results than noise- 
assisted, essentially classical transport processes, though 
requires additional optimization of the molecular struc- 
ture. Indeed, this even holds in the presence of not too 
strong ambient noise, and thus identifies yet another sce- 
nario where genuine quantum effects may define an evolu- 
tionary advantage. 

Random molecular networks. The noise- 

enhanced as well as the multi-path-interference scenario 
start out from the same structural elements that are un- 
ambiguously given by experimental observations: Exciton 
transfer occurs across a random molecular network with 
local sites effectively modeled as electronic two-level sys- 
tems, and further background degrees of freedom, which 
may possibly exert some effective driving 
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principle be accounted for by some non-Markovian envi- 
ronment coupling 13 . Experimentally observed exciton 
cross-terms [l] , though with quite some experimental scat- 
ter [2] , suggest excitonic coherence times which exceed the 
exciton transfer time between antenna and reaction center 
by a factor two to five [l], even at ambient temperatures 
[2]. The presently best available effective Hamiltonian, 
which is inferred from experimental data and advanced 
quantum chemical model calculations predicts a strong 
suppression of strictly coherent transport across the FMO 
complex (approx. 5 % transfer efficiency 2Tp2] ), from the 
input to the output site, due to predominantly destructive 
multi-path interference upon transmission. 

However, exciton transfer in the FMO complex occurs 
with efficiencies larger than 95 % [t], and there are es- 
sentially two alternative scenarios which can explain this 
discrepancy between the experimental data and the best 
microscopic Hamiltonian presently available. Since quan- 
tum coherence can be destroyed by ambient noise, and 
since suppressed transport under purely coherent dynam- 
ics can only be due to destructive interference effects, it is 
very natural to argue in favor of noise-induced transport 
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Since, however, destructive multi-path interference 
is known to be very sensitive with respect to changes of 
the boundary conditions and/or the Hamiltonian's cou- 
pling matrix elements [24| , and since even the best avail- 
able model Hamiltonian for the FMO complex is garnished 
with appreciable uncertainties for its individual entries 
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one may equally well argue in favor of constructive 
multi-path interference as the observed efficiency's cause. 
This even more so since all experimental data that so far 
lend support for the coherent transport hypothesis are ob- 
tained from bulk measurements rather than from single 
molecule spectroscopy, and therefore might mask much 
longer coherence times by inhomogeneous broadening ef- 
fects [3]. 



Since both alternative explanations allow to predict 
large transfer efficiencies, let us have a closer look at the 
respective key mechanisms. We model the energy con- 
serving, unitary dynamics of a single excitation injected 
into a molecular network alike the FMO complex by the 
Hamiltonian 



H 
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(1) 



with \i) and \j) the electronic states where the excitation 
is localized at the individual molecular sites i and j, re- 
spectively. We assume that initially only one site, "in" , is 
excited, which is identified with the first site of the net- 
work. The sites No. 2 to 6 are referred to as the "interme- 
diates" . The seventh and last site is the designated output 
site, "out" , where we add an energy sink, such as to mimic 
the irreversible exciton absorption at the reaction center. 
We couple each of the seven molecular sites to a private 
(i.e., there is no inter-site communication through the en- 
vironment) dephasing environment. Sink and dephasing 
induce some irreversibility on the FMO degrees of free- 
dom, which we incorporate by the Lindblad terms 

isi„k(e) =r(|0)(out|e|out)(0| - i{|out)(outU}), (2) 

where |0) and {, } are the ground state of the molecular 
network and an anticommutator, respectively, and 



(3) 



into the effective evolution equation of the excitonic state 
on the network, 

g{t) = -i [H, g{t)] + Lsink(e(t)) + L^cpiMt)) (4) 

(where we have set h = 1 for convenience) . To obtain a ro- 
bust comparison of the different transport mechanisms, we 
statistically sample the transport efficiency over different 
realizations of H, by random sampling over the positions 
of all intermediate sites, within a sphere with input and 
output site placed on its north and south pole, respec- 
tively. Random positions Vj translate into random real- 
izations of H by defining the coupling matrix elements as 
a function of the intersite distances. 



(5) 

I . By choosing ran- 



with some constant a, and r, 
dom positions for the intermediate sites, we in some way 
take into account the experimental uncertainty concerning 
the actual FMO Hamiltonian 
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Let us stress, however, 
that we do not present our model as a realistic description 
of a particular experiment. On the contrary, we are rather 
interested in general properties of transport, as they will 
appear if we allow for a large variety of possible configu- 
rations, without imposing, from the very beginning, strict 
experimental boundary conditions. 
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The model discussed above has three intrinsic time 
scales which will largely determine the expected transfer 
efficiencies — the direct exciton transfer time 



T = 



(6) 



between input and output site, in the absence of all inter- 
mediate sites, the local dephasing rate 47 (identical for all 
sites), and the sink dissipation rate F, which we will fix at 
the value F = 10/T in the following Randomly placed 
intermediate sites between input and output help to en- 
hance or suppress coherent transport (21 22 , by larger 



coupling strengths between closer sites and appropriate 
phase relationships. The dephasing rate defines the time 
scale Tdeph — (47)^^ on which such phase relationships 
can have a bearing on the overall transport behavior, and 
the sink dissipation rate defines an optimal time scale on 
which population has to be delivered to the output site, 
to make it immediately available for the reaction center. 

Statistics of transport efRciency. To compare 
the transport efRciency provided by different molecular 
conformations and different dephasing rates, we define the 
average transfer time 



tPoutit) dt, 



(7) 



which is the expectation value of the time required to 
transfer the excitation to the sink, determined by the pop- 
ulation Pout(i) = (0|^{i)|0)/F of the output site. We see 
here that the ground state |0) can only be populated by 
delivering the exciton from the output site to the sink, 
with rate F, cf. also eq. ([2]). The shorter the transfer 
time, the more efficient the transport 

We are now prepared for a statistical analysis of exci- 
ton transfer times across a molecular network alike the 
FMO complex, to assess the potential role of coherent vs. 
noise-assisted transport mechanisms to steer transport ef- 
ficiencies. In order to draw a landscape of the exciton 
transport efficiency in molecular networks as modeled by 
eq. Q, we sample over fifty million random and distinct 
conformations. For each conformation in the ensemble 
the computational procedure involves the following tasks: 
First, the positions of the five intermediate molecular sites 
are randomized, which are then used to populate the ma- 
trix entries of the Hamiltonian H according to eq. ([s]). 
For values of 7 from 10~^F to 10"^r in 200 logarithmic 
steps, i.e. for essentially vanishing to very strong dephas- 
ing rates, the master equation, eq. Q, is solved via exact 

^This implies one incoherent event on the time scale T/10, hence 
a rather efficient drain. We have used this time scale as an efliciency 
benchmark in earlier work [21[|22| . 

■^Note that X is a reasonable efliciency quantifier if efficiency is 
qualified as rapid and irreversible excitation transfer to the sink, 
but that it does not distinguish quantum from classical transport 
efficiencies, since it integrates over all times. It is however evident 
that the definition of any efficiency quantifier is a matter of prag- 
matic choice rather than of principle, and that all such quantifiers 
call for a careful interpretation 



numerical diagonalization and the transfer times T are cal- 
culated. The last step employs logarithmic data binning 
to record the transfer time histogram fi, shown in Fig.[l| 
as a function of 7. In the left density plot of the figure, 
the grey scale represents the probability density of con- 
figurations giving rise to a certain transfer time, under a 
given dephasing rate. Configurations above the left hand, 
decreasing dotted line correspond to transfer times longer 
than Tdoph, j-e. the exciton transfer is here due to classical 
transport across the network. For even larger dephasing 
rates, transport is "frozen" due to a Zeno-like projection 
mechanism (Tzono oc 7, as indicated by the increasing dot- 
ted line), while configurations below the left hand dotted 
line achieve efficient transfer on time scales shorter than 
"Jdcph- The latter thus mediate coherent exciton transfer 
and are the only ones which are eligible for claiming an 
unambiguous quantum enhancement of exciton transfer in 
the FMO complex. Indeed, the shortest transfer times 
are observed for a finite subset of these configurations, on 
the lower left hand side of Fig. [T] If we furthermore op- 
timize the molecular configuration in the absence of any 
dephasing, for the same sink dissipation, and then expose 
this optimal configuration to finite dephasing rates, we 
obtain the dashed curve in Fig. [T] which exhibits exciton 
transfer times shorter than the classical value, associated 
with the highest density of configurations in the plot, by 
at least a factor two, for almost all dephasing rates that 
allow for coherent transfer. Note that this optimal config- 
uration still yields very efficient excitation transfer even 
for higher dephasing rates, which suggests that the opti- 
mized coherent transport on short time scales imparts an 
initial advantage, even when classical activation takes over 
on long time scales. 

Discussion. — The optimization landscape thus pro- 
vides a clear picture of the possible strategies to achieve 
efficient transfer across an FMO like, disordered network: 
The large majority of randomly sampled configurations 
requires the assistance of noise to allow for efficient trans- 
port, which, however, will be classical in nature, and will 
occur on time scales clearly longer than Tdeph- Indeed, 
the clustering of all the configurations associated with this 
transport mechanism around T ~ 10/F = T and 7 ~ F 
in Fig. [1] expresses an approximate matching condition of 
dephasing and sink dissipation rates, which is an intrin- 
sically classical, statistical synchronization phenomenon. 
Excitation transfer here stems from the noise induced de- 
struction of quantum coherence, is induced on time scales 
longer than the system inherent coherence times, and is 
essentially independent of the microscopic Hamiltonian 
which generates the coherent dynamics on the network 
(hence the large density of configurations at these trans- 
fer times and dephasing rates). Under an evolutionary 
perspective, this non-selectivity with respect to the mi- 
croscopic coupling structure can therefore neither define 
an evolutionary advantage by optimizing the molecular 
structure — most structures will do, as also highlighted by 
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Fig. 1: (a) Probability density fi of the average excitation transfer time 1, eq. ([7|, for N — 7 molecular sites and sink rate 
r — lO/T, as a function of the dephasing rate 7. The two dotted, diagonal lines are given by the dephasing time Tdoph = (47)~^, 
and by an approximate Zeno time Tzeno oc 7, respectively. On time scales T > Tdeph, the purity of the excitonic state on the 
molecular network has dropped to its minimum value, hence the transport is essentially classical. The white line shows the 
median T, the dot-dashed line the minimum transfer time, and the dashed line the transfer time of a configuration that has 
been optimized for 7 = 0. (b) Height profiles of the probability densities for fixed dephasing rates 7 = (solid line), 10~^r 
(dashed line), 1.06r (dotted line), lO^F (dash-dotted hne). 
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the pronounced minimum of the median of the distribu- 
tion (white curve) at 7 ~ F. 

In contrast, Fig. [l] also shows that a smaU but finite 
sub-ensemble of optimal configurations mediates efficient 
quantum transport of the excitation from input to output, 
faster than the noise-assisted, classical transport, and due 
to constructive quantum interference upon transmission, 
which prevails in the presence of noise, actually in the 
entire range of dephasing rates considered in Fig. [ij If 
quantum mechanics is at the origin of the experimentally 
observed exciton transfer efficiencies in the FMO complex, 
it therefore must stem from such type of optimal molec- 
ular configurations, or otherwise ought to be induced by 
non-Markovian 14 or driving effects 18 . Since the opti- 



figurations achieving maximum performance due to con- 
structive quantum interference will certainly spur the de- 
sign of new experiments and, in the long run, advanced 
devices such as a new generation of organic solar cells that 
utilize the beneficial aspects of quantum coherence. 
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